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Abstract. Fractional control systems with the Caputo derivative are consid- 
ered. The modified controllability Gramian and the minimum energy optimal 
control problem are investigated. Construction of minimizing steering controls 
for the modified energy functional are proposed. 



1. Introduction 

The Fractional Calculus is an important Mathematical discipline [TBI HH1 HI]. 
Several recent books on the subject have been written, illustrating the usefulness 
of the theory in applications [16j 1201 [211 [26] ■ Different notions of fractional- 
order derivatives are available, including the Riemann-Liouville, the Griinwald- 
Letnikov, and Caputo, as well as the generalized functions approach. One of the 
youngest fractional derivative, formulated in 1967, is the Caputo derivative. The 
main advantage in using the Caputo derivative is that it avoids problems connected 
with initial conditions of fractional differential equations. Indeed, when an initial 
value problem is formulated for the Caputo derivative, initial conditions look like 
in the classical (non-fractional) case. 

The state-space description of fractional order control systems is developed 
in [71 [H [H OH [21 [H]. Controllability and observability of finite-dimensional 
fractional differential systems are investigated mainly in |17j . for systems with 
Riemann-Liouville derivatives. For positive continuous-time linear systems, the 
reachability property has been worked out in [14] . Differently, here we deal with 
the controllability property of finite-dimensional linear fractional-order differential 
systems via the Caputo derivative. The novelty is the construction of control laws 
corresponding to the minimal modified energy, where we use a neutralizer of the 
singularity similarly to the one defined in [17] , 

The paper is organized as follows. In Section [2] we recall the main properties of 
the generalized Mittag-LefHer function and the fractional analogue for the exponen- 
tial matrix. We also review some properties of fractional integrals and derivatives. 
In Section [3] we discuss the Gramian formula for fractional control systems and we 
prove a new result about the minimal modified energy steering control law (The- 
orem [O]). Illustrative examples are discussed. In Section [4] we recall the classical 
conditions for controllability and propose the construction of a different control law 
for fractional control systems (Theorem 14.31) . 
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2. Preliminaries 

We investigate control systems with commensurate order a G (0, 1]. For a = 1 
the classical situation is obtained; the novelty appears when a € (0, 1). 

Definition 2.1 ([IU]). We denote by E a ^ the two-parameter Mittag-Leffler func- 
tion defined by the series expansion 

OO 

(2-1) ^W = Er7TTm' a>0 ' f3>0 - 

When ft = 1, we use the notation E a = E a<1 . Let A e K" xn . We extend (|2~T|) to 
the matrix case as follows: 



, „ r(fea + /3) 
Definition 2.2 (QS]). Let A € M nx ". By 

e At = ta -i gA* r[(fc + 1)a] = g^ r[ (fc + l)a] = ^-(A") 

we denote the a-exponential matrix function. 

For a = 1 we have £i(At) = e^* = exp(^4t), where exp denotes the classical 
exponential matrix. We remark that the properties exp(A + B)t = exp(At) exp(Bt) 
and exp~ 1 (At) = exp(— At), satisfied for square matrices A and B, are not valid 
for functions e^* and E a (At a ). 

Definition 2.3 ([M EH HQ). Let tp e Li ([t ,ti],M). The integrals 



r(«) 



ig_V»(t) = t^t / 1 V(r)(r - t^^dr , t < t x , 

r(a) y t 

where T is the gamma function and a > 0, are called, respectively, the left-sided and 
the right-sided fractional integrals of order a. Additionally, we define the identity 
operator I by / := 7 t ° + = J t °_. 

We have the following integration by parts formula for fractional integrals. 

Proposition 2.1 (16;). Let a > and 1/p + 1/q < 1 + a, p > 1, g > 1, 

p ^ 1 and i;^ 1 in t/ie case 1/p + 1/g = 1 + a. Then, for tp € L p ([to,ti],R) and 
(/> G L 9 ([to,ti],R), the following equality holds: 

(2.2) r ^(r)/r o+ ^(r)dr = T ^/^(t^t . 

Jt "'to 

Definition 2.4 ([HE US])- Let </? be defined on the interval [to,ti]- The left-sided 
Riemann-Liouville derivative of order a with lower limit to is defined by 

< 2 ' 3 ' D '°+*M " whj (s)"X'" ( ™' - tT "" Ut ' 

where n is the natural number satisfying n = [a] + 1 ( [a] denotes the integer part of 
a). Similarly, the right-sided Riemann-Liouville derivative of order a with upper 
limit ti is defined by 

1 / r\ \ n Z"* 1 

(2.4) ^) = f(^)(-^ I <P(T)(T-tr-°-^. 



Minimal modified energy control for fractional linear control systems with the Caputo derivative 3 



The next proposition is based on [3TJ Corollary 2, p. 46] and is particularly useful 
for our purposes. 

Proposition 2.2 (Integration by parts). Let f g I^_(L p ) and g € I" g+ (L q ) with 
l/p + l/q < I + a. The following formula holds: 



(2.5) / f(t)D? o+ g(t)dt = / g(t)Dl_f(t)dt, < a < 1 . 

Jt Jto 

Proof. Denote D*_f(t) = ip(t) and D« a+ g{t) = ip(t). The equality follows 
from (E3 since I? 0+ D? o+ f(t) = f(t) is valid for / e (cf. fS]). □ 

In the description of fractional control systems we use the notion of Caputo 
derivative, which is the preferred fractional derivative among Engineers. 

Definition 2.5 ([IS]). The left- and right-sided Caputo fractional derivatives of 
order a > on [to,ii], denoted by c 'Df g+ ip(t) and c D" i _ip(t), respectively, are 
defined via the Riemann-Liouville fractional derivatives (|2.3[) and (|2.4p by 



and 

(,T) °*M t ) = -*,-(rt)-f ^ lt '*'-* ), 

where n = [a] + 1 for a £ No and n = a for a € No- 

When < a < 1 the relations (|2.6p and (|2.7p take the following form: 
C D? Q+V (t) = D? o+ (<p(t) - <p(to)) , C DlMt) = K- (*>(*) ~ • 

Let AC[io,ii] be the space of functions that are absolutely continuous on [tojii] 
and AC n [to, t{\ denote the space of functions ip that have continuous derivatives up 
to order n — 1 on [to, t\] and such that ip 1 - 71 ^ 1 ^ g AC[a, b\. 

Proposition 2.3 ([IB]). Let a > 7 n = [a] + 1 j/a N and n = a i/aeNj. J/ 
iy2 6 AC"[to) ii]) t/jen t/ie Caputo fractional derivatives c Df + tp(t) and c D" i _ip(t) 
exist almost everywhere on [in,ii]. Moreover, 
(a) If a ^ N , then 

C DZM*) = f^Taj [ V (n) (T)(t - rr—'dr = J*- (*,<»>) (t) 



to 



f&J Ifa = nG N , then c D? 0+ ip(t) = ip (n ^(t) and c D^_(p(t) = {-l) n (p^(t). 
For a function a; : [0, T] — > M" we use similar notation as in the classical case: 

/ Xl (t) \ ( C D% +Xl (t) 



C D% + x{t) = C D^ 

\ x n {t) ) \ c D^ + x n {t) 



Such situation, when for each component we use the same fractional order a of dif- 
ferentiation (in the Riemann-Liouville or Caputo sense), is known in the literature 
as the fractional derivative with commensurate order [14j [22] . 

Proposition 2.4. For a > the following holds: 
(i) c D? o+ E a (A{t - to) a ) = AE a (A(t - t ) Q ); 
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(ii) Df 0+ e^-^=Ae^-^; 

(lii) D°_S(T - r) = AS(T - r), where S(t) = ef. 

Proof, (i) Directly from the definition of the classical Mittag-Lemer function E a , 
and from the formula of the Caputo derivative of a power function, one has: 

Since the Caputo derivative of a constant function is zero, 

v (+ 4- \ka v (+ + \(k-l)a 

= A£ Q (A(t-t ) a ). 

(ii) As the formula for the a-exponential function does not consist of constant 
terms, and since the formula for the Riemann-Liouville derivative of a power func- 
tion is the same as for the Caputo derivative, we have: 

where we used the fact that lim jJ-^ = 0. 

(iii) Using the formulas 

(T - rf- 1 = (T - r)^"- 1 , lim D%_ (T " ^ - " 



r(/?-a) v ' ' /3^« J " r(a) 



(cf., e.g., [21]), we get: 



^5(T - r) = D ? _ (/^(T - .)- + 

= AP , -;r i + ^ (r 7 ) r+- 

r(a) r(2a) 

1 j ^ r[(fc + l)a] 
= AS(T - t) . 



Proposition 2.5. For a > the following relation holds: 

E a (A(t - t ) a ) = I + ( Aei^dT. 

Jt 

Proof. Follows by direct calculation of the integral: 

ft ft °o \(A+l) a -l 

=±^^- = Ea{ A (t - to r)-i. 

k—l ' 



□ 



□ 



Since e At = t a 1 E a ^ a (At a ) and each Mittag-Lemer function E a _ a (az a ), a > 0, 
is an entire function on the complex plane, we can state the following: 

Proposition 2.6. Let a > 0. There is a uniquely determined function g(t) = 
t 1 - Q G(i) such that e At g(t) = E a , a (At a )G{t) = I, for t ^ 0, and \\me At g{t) = I. 
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Lemma 2.1. Let a > and i/j(t) £ lQ + (L q ), where l/p + 1/q < 1 + a for p such 
that all components of S(T — t) belong to I^_(L p ). Then, 

S{T-T)D« + lP(T)dT = ( AS{T - T)4>{ T )dT . 



>o Jo 
Proof. Taking into account (I2.5[) and item (iii) of Proposition 12.41 we have 

S(T - T)D% + ij){T)dT = [ ^(t)D t _S{T - r)dr = f AS(T - T)^(r)dT . 



a 



3. Control systems and the Gramian 



We consider the following linear time-invariant control system of order a £ (0, 1], 
denoted by S: 

c Dg + x(t) = Ax(t) + Bu{t) , y = Cx(t) , 

where x(t) £ R™, u(t) G R m , matrix A € M nxn (R), B € M nxm (R), C £ M pxn (R), 
and Dq + indicates the fractional Caputo derivative of commensurate order a. 

The forward trajectory of the system S, starting at to = and evaluated at 
t > 0, is the solution of the initial value problem c DQ + x(t) = Ax(t) + Bu(t), 
x(0) =o6l" [IS]: 

(3.8) 7(t, a, u) = f I + ^ S*(r)AdT^ a + ^ S(t - t )Bu(t )dr , 

where S'(i) = e^*. Moreover, we can represent (|3.8[) in the following way: 

7(t, a, u) = So(t)a + / S(t — T)Bu(r)dT , 
Jo 

where Sa(t) = -E a (At Q ) = / + J * 5(r)A(iT. The formula for the forward trajectory 
can be obtained using the Laplace transform [14.. Taking into account the output 
of E, the forward output trajectory is then defined by values evaluated at t > 0: 

r?(t, a, u) = C-y(t, a, u) = C (i + J S(T)Adr^J a + C S(t- t)Bu(t)(1t . 

For system £ we define the notion of controllability in the standard manner: 

Definition 3.6. Let T > 0. The system £ is controllable on [0, T] if for any a £ R" 
and b £ 1" there is a control u(-) defined on [0, T] which steers the initial state 
7(0, a,u) — a to the final state j(T, a, u) = b. 

Following [T71 [25] we denote by 

Qt= S(T -t)BB*S*(T -t)(T -t) 2{1 - a) dt 
Jo 

the controllability Gramian of fractional order a, on the time interval [0, T], cor- 
responding to the system E. As in the classical case [35J, Qt is symmetric and 
nonnegative definite. The term (T — t) 2 ( 1 ^ Q ) under the integral is called in [17] a 
neutralizer of the singularity at t — T. It is needed in order to ensure the conver- 
gence of the integral. 

Let T > 0. By L 2 a ([0,T],R m ) we denote the set of functions ip : [0,T] -> R m 
such that tp defined by tp(t) = (T — t) a ~ 1 tp{t) is square integrable on [0, T]. 

Theorem 3.1. Let T > and Qt be nonsingular. Then, 
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(a) for any states a, b € W 1 the control law 

(3.9) u{t) = -(T-t) 2{l -^B*S*{T-t)Q T 1 f T (a,b), te[0,T), 

where 

f T (a, b)= ^1 + S(t)Adtj a-b=-b+ S {T)a 

and u{T) = 0, drives point a to point b in time T; 

(b) among all possible controls from L\ ([0,T],K m ) driving a to b in time T, 
the control u defined by \3. 9\) minimizes the integral 



(3.10) ( T \{T ~t) a - l u{t)\ 2 dt. 

Jo 

Moreover, 

T 

|(T - tr^uit^dt =< g T Vr(a, b)J T (a, b) > , 

where < • , • > denotes the inner product. 
Proof, (a) From the expression of u we have that 
j(T,a,u) = fr{a,b) + b 



(jf (T - t) 2 ^- a) S(T - t)BB*S*(T - t)dtj Q T 1 / T (a, b) 



= f T (a, b) + b- QtQt 1 !^ b) = b. 
Moreover, lim u(t) = u{T) since lim (T - t) 2{1 - a ^ B* S*(T ~ i)Q T x f T (a,b). 
(b) We notice that for h(t) = \(T - t)\ 2 ^~^ 

/ \{T -t) a ~ 1 u{t)\ 2 dt= / \{T-t) 1 - a B*S*{T~t)Q^ 1 f T (a,b)\ 2 dt 
Jo Jo 

r T 

= / h(t)(B*S*(T -t)Q^ 1 f T (a,b),B*S*(T ~t)Q^f T (a,b))dt 
Jo 

= ^jf h(t)S(T - t)BB*S*(T - t)dt, Q^f T (a, b)\ 

= (QtQt 1 / T (a ! 6),g T 1 Ma,b)) = (f T (a,b),Q^f T (a,b)) . 

Let us take another control u for which (T — t) a ~ 1 u(t) is square integrable on [0, T] 
and 7(T, a, u) = b. Then, 

T 

( T _ t )2(a-1) ( u ( t ) t u( t ))dt 



Hence, 



ft(t) (u(t), (T - t) 2i - 1 -^B*S*{T - t)Q^f T (a, b)\ dt 

[ (u(t), B*S*(T - t)Q^f T (a, b)) dt = (f T (a, b),Q T l f T {a, b)) . 
Jo 

(T — i) 2 ( 1_Q ) (u{t),u(t)) dt = / (T -i) 2 ^-") (u(t),u(t)) dt 
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and from that we obtain 

{T-tf^-^\u{t)\ 2 dt 

o 

,T ,T 

= / (T-t) 2 ^-^\u(t)\ 2 dt+ / (T - t) 2(1 ~ a ^\u(t) -u{t)\ 2 dt , 
Jo JO 

which gives the minimality property for the integral. □ 
As in the classical case, the operator 

u^j{T,0,u)=[ S(t)Bu(T - t)dt 



is a linear operator from the space L 2 a ([0,T],R m ) into R n . Hence the classical 
result can be replaced by the following proposition. 

Proposition 3.7 (cf. // any state b £ M. n is attainable from a — 0, then the 

matrix Qt is nonsingular for any arbitrary T > 0. 

Example 3.1. Let £ be the following system evaluated on M 2 : 

( c D°4 Xl (t)=x 2 (t), 
• \ c D^x 2 (t)=u(t). 

Let us take a = (1,0)* and b = (0,0)*. Since A = ( ° Q * ^ and B = f ® 

we obtain that S(t) = ^ ^ J_ ^ while S (i) = ^ q ^- Hence ' the 

formula for the solution with the initial condition 7(0, a, u) = a is 

7 (M,»)=( J f )« + /'( ^) B " M<fr ' 

Let us take w(t) = 1. Then, for the given a, ^(t, a,u) = ^1 + t 2-^=Q . From the 

last expression we see that using constant u(-) = 1 for t > we are not able to steer 
the given initial point a to the origin. 

Let now /r(a, b) = Sa(T)a — b = a. The Gramian has the form 



2 3^ 



and the control 



y2 1 y3/2 

drives a to b with the modified energy 
m 







\(T-t)-^uit)\ 2 dt = ^ 



Example 3.2. Let a £ (0, 1). Consider the following fractional system £ on M 2 : 

| c ^ + x 2 (i) = - Xl (t)+u(t). 
The matrix A is now skew-symmetric, and thus 

.1" /• S i _° J I .!r 
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Hence, A k = I if k = 0,4, 8, . . .; A k = A if k = 1,5,9,...; A k = -I if k 
2, 6, 10, . . .; and A k = -A if k = 3, 7, 11, . . .. Moreover, 

S(t) = t"- 1 J-±- + A— j—- - 1^-— - A- 



T(a) r(2a) r(3a) F(4a) 

j-Q — 1 ^-3a — 1 \ / y.2a— 1 jAch—I 



v T(a) r(3a) + J + A V r (2") r ( 4 ") 



Using the notation 



£2(fc+l)a-l 00 ^ ^(2fe+11o:-l 



sin Q t = V(-l) fe -__ _ -, cos Q i = V(-l) fc - 



^ r[2(fc+l)a]' ^ ' r[(2fc + l)a]' 

„,,, / cos Q t sxn a t \ . _ / \ , 
we can write Mi = . , , . As B = , , we nave 

y — sm„ t cos Q t J \ 1 J 

Qt= [ (T - t) 2(1 - Q) • M a (t)dt 



with 

M a (t) = 



sm 2 a (T-t) sin Q (T - t) cos a (T - t) 

sm a (T - t) cos Q (T - t) cos 2 a (T~t) 



To get an exact formula for Qt and Q^ 1 is difficult. We can, however, easily obtain 
approximations with the desired precision for both matrices, and an approximate 
formula for the optimal control. Let us consider a concrete situation. Let a = 5. 



1 / +2? 2 k e k - x \ 1 , 

cosi t = -= 1 - > —7 = -=(1 -2i + 

We choose to approximate cosi t by functions 

c L {t) = ^=[l-Y j ^ , L = l,2,.. 

In this way we can approximate Qt by 

-2(T-i) , T _.s (r-t) 



[ j( Ci (T-t)e-( T -*) cf(i) ' 



Let T = 10, the starting point to be a — (0 1)*, and the final point to be the 
origin, i.e., b = (0 0)*. Applying formula (|3.9[) for the steering control u with 
the approximate expressions for S*(T — t) and Q^. , we check that we are very 
close to the final goal b. Taking more terms in the functions cl, the value of 
the modified energy tends to the minimum. For example, for L = 1 we have 
c\(t) = "^t=|(1 — 2t) and the approximate value of energy is mi = 1.02; for L = 11 
we have ran — 0.0921; and for L — 12 we get m\2 = 0.0911. These numbers were 
calculated using MATLAB with a final precision of four digits. 



4. Rank conditions and steering controls 

In the classical theory it is quite easy to explain why the condition rankB = n 
is sufficient for controllability. Roughly speaking, one needs to construct a special 
control u(-). In the case of fractional order systems we use the special function 
given by Proposition ^. 61 
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Proposition 4.8. Let rankB = n and B + be such that BB + = I. Let g(-) be the 
matrix function defined by Proposition \2.b\ Then the control 

u(t) = ^B+g(T -t)(b- S {T)a) , t G [0, T] , 

transfers a to b in time T > 0. 

Proof. The proof follows by a direct calculation: 

i r T 

7(T, a, u) = S (T)a +- S(T - t)BB+g(T - t) (b - S {T)a) dt = b . 
1 Jo 

□ 

Example 4.3. Let E be the following fractional system evaluated on R: 

E : C D* + x(t) = u(t) , 
where a E (0, 1]. Let us take a, b g R, T > 0. Then, 



S(t) = —, S (t) = l, B = B^ 
r(a) 



and 



T 



7(T, a, u) = a + — L / (T - t) a - l u(t)dt . 
r(a) 



(i 



Consider, accordingly to Proposition 12. 61 <?(£) = Q T(a). From Proposition 

u(t) = ^(T-t) 1 ~ a (b-a) 

transfers a to 6: j(T, a, u) = b with the modified energy p.lOj) of value 

T 2 (a)(b-a) 2 

m = . 

T 

It is also the minimum energy as u{t) — u(t). 

An algebraic condition equivalent to controllability for fractional linear control 
systems has been derived in [T3] and cited again in [7J[T7]- According to Theorem l3.ll 
and Proposition 13.71 we can state a similar formulation as in |25] : 

Theorem 4.2. The following conditions are equivalent: 

(a) An arbitrary state b € R n is attainable from 0. 

(b) System E is controllable. 

(c) System E is controllable at a given time T > 0. 

(d) Matrix Qt is nonsingular for some T > 0. 

(e) Matrix Qt is nonsingular for an arbitrary T > 0. 

(f) rank [A\B] = rank [B, AB, . . . , A n ~ 1 B~\ = n. 

If the rank condition is satisfied, then the control u(-) given by (|3.9p steers a to 
b at time T. Our goal now is to find another formula for the steering control by 
using the matrix [A\B] instead of the controllability matrix Qt- It is a classical 
result that if rank [A\B] = n, then there exists a matrix K € M(mn,n) such 
that [A\B] K = I g M(n,n) or, equivalently, there are matrices K-y, K2, ■ . ., 
if„ G M(m, n) such that Bifi + ABK 2 H h A n ~ x BK n = L. 

For the next construction it is more convenient to use the notion of Riemann- 
Liouville derivative. We begin by introducing a notation for compositions of the 
Riemann-Liouville derivatives with the same order a, a G (0,1). Let i?o+W) = 

ip(t). Then for j G N, recursively, we put i?o+ + V(*) := ^0+ (#o+V>(*)} 
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Theorem 4.3. Let rank [A\B] — n and a £ (0, 1). Let p be such that S(T — t) £ 
It—(L p ) and tp be a real function given on [0,T] such that 

(i) fiip{t)dt=l; 

(ii) Rfim G Io+( L i) forj = 0,..., n-1, where ${t) = g{t) (b - S (T)a) <p(t) 
and S(T-t)g{t) = I,t£ [0, T], for 1/p + 1/q < 1 + a. 

Then the control u(t) = K^t) + K 2 D% + tp(t) H h -K'n-Ro-f^VW, t £ [0,T], 

transfers a to b at time T > 0. 

Proof. Using j — 1 times the formula (|2.5p of integration by parts and Lemma 12.11 

f S{T-t)BK j R^ 1 ij)(t)dt= j SiT-tfA^BKj^dt 
Jo Jo 

for j = 1, . . . , n. Then, 

S{T-t)Bu(t)dt= [ S(T - t)Bip{t)dt 
Jo 

and 7 (T, a, u) = S (T)a + ^ S{T - t)g{t) (b - S (T)a) ip(t)dt = b. □ 

5. Conclusion 

The concept of fractional (i.e., non- integer) derivative and integral is increasingly 
being recognized as a good tool to model the behavior of complex systems in various 
fields of science and engineering [H 131 SI El HI HOI H3] ■ This is particularly true in the 
areas of control theory and control engineering, where fractional order controllers 
and plants provide confirmed evidence of better performances when compared with 
the best integer order controllers [H El El El [HI H21 [HI [13 HZ]- In this work we 
consider linear time-invariant control systems of fractional order a, a £ (0, 1), in 
the sense of Caputo. The controllability question for such systems is addressed. 
Main results give sufficient conditions assuring the existence of controls that drive 
the system between any two desired states. Explicit formulas for such controls 
are given. A control law is obtained that minimizes the modified energy J Q T |(T — 
t) a ~ 1 u(t)\ 2 dt. Our results are reduced to the classical ones when the fractional 
order of differentiation a tends to one. 
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